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The next-to-leading order QCD corrections to the Standard Model process 

"X' 

of single top production, bg tW , for the CERN Large Hadron Collider with 
y/s = 14 TeV have been calculated. For renormalization and factorization scales 
(jl = no (fj,o — m t + mw)i the NLO hadronic cross section is ~ 37 pb, while ~ 25 
pb for tree level. The NLO QCD corrections can enhance the cross section by a 
factor from 1.33 to 1.66 for ^ < \i < 2^o- 
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The top quark plays an important role for testing the standard model (SM) and search- 
ing for new physics beyond the SM, due to its large mass, the same order as the electroweak 
symmetry breaking scale. In order to carefully measure the top-quark electroweak interac- 
tions it is useful to consider single top production, in addition to studying the decay of the 
top quark in tt events. Within the context of the SM, single top production modes provide 
a direct measurement of the Cabbibo-Kobayashi-Maskawa matrix element Vtb- 

At hadron colliders, single top quarks can be produced within the SM in three different 
channels, the s-channel W* production Jl|-|5|, the t-channel W-exchange mode fl5Hl!^|, and 



through tW~ production [|13],|rjJ. These three subprocesses have very different kinematics 
and experimental signatures, and are sensitive to different types of new physics in the top 
quark sector [[nj. It should be noticed that the tW~ production rate is extremely small at 
Tevatron, but is much greater at LHC The study shows |T6| that a 5cr observation of tW~ 
signal is possible at very low luminosity at LHC, with 20 fb^ 1 , cross section can be measured 
up to the accuracy of 1%. 

Both the precise measurement of Vtb arid the indirect detection of new physics require an 
accurate calculation of the single top quark production cross section. The QCD corrections 
to the s-channel W* production ]3| and the t-channel W-exchange mode |ni|I5[ have been 
done. However, up to now, only part of the QCD corrections [0(1/ log mf/mf)} to the cross 
section for pp -+bg — *> tW~ has been known | [ 115 , 171 . In this letter, the results of the 



complete next-to-leading-order QCD correction to tW~ production will be presented. A 
detailed review of the calculation will be published elsewhere Ul8| . 

The Feynman diagrams for tW~ production via the parton process b(pi)g(p2) — > 
t(k\)W~(k2), including the QCD corrections, are shown in Fig. 1. The Born diagrams are 
shown in Fig. |l](a), the NLO diagrams by virtual gluon-exchange and the gluon-radiation 
(gr) are shown in Fig. [l|(b) and Fig. |l](c). Because the topologies of the initial-gluon (ig) 
diagrams are the same with the gluon-radiation processes, we don't show here the diagrams, 
which can be easily obtained from Fig. |TJ(c) by treating two gluons as initial partons. The di- 



tln Ref. [17|, the QCD corrections to the similar process of Wc production at Tevatron has been 



calculated. 
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agrams are created by use of FeynArts fl9l , |2~0|j and are handled with the help of FeynCalc |21 . 



We perform all the calculations in d = 4 — 2e dimensions and adopt MS renormalization 
and factorization schemes. 

The total cross section for pp — > bg — > tW~ at 0(a 2 ) can be written as: 

a( s ) = a\s) + a vir {s) + a 9r {s) + a i9 (s), 

= a 2body (s) + a 9r (s) + a l9 (s), (1) 



where 



a 2 body (s) 



dL\ f 1 dx 



,z 



2 



\dz) b = 2z ] z2 ^r fb / p ( x ^f)f9/p(—,H) ( 2 ) 

with zq = (mw + mt)/\/s. Here cr°, a mr , cr 9T ' and a* 9 are contributions from tree level, 
virtual, gluon-radiation and initial-gluon diagrams. The two-body subprocess cross section 
can be expressed as 

a 2bo dy{§)= jY^\ Mren \i d <s> 2 

= jY}M Q \ 2 d<& 2 + jY,2Re{M vir M+)d<t> 2 

= a° + a mr . (3) 

Here Mo, M ren and d$2 are the tree level amplitude, the renormalized amplitude and the 
two-body phase space in d dimension. The details of the renormalization procedure and 
the explicit expressions of M ren will be given in Ref. |Tj|. As usual, a mr contains infrared 
divergences after renormalization, which can only be cancelled by adding contributions from 
a 9r remaining collinear divergences are absorbed by the redefinition of the parton 

distribution functions (PDF). 

The real corrections a gr and a %9 have been computed using the two cut-off phase space 
slicing method (TCPSSM) |[22|| . The main idea of TCPSSM is to introduce two small 
constants 5 S , 5 C . The three-body phase space can then be divided into soft and hard regions 
according to parameter S s , and the hard region is further divided into collinear and non- 
collinear regions according to parameter 5 C . In the soft and collinear regions, approximations 
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can be made and analytical results can be easily obtained. In the non-collinear region, 
numerical results can be calculated in four dimension by standard Monte Carlo packages 
because it contains no divergences. The physical results should be independent on these 
artificial parameters S s and S c , which offers a crucial way to check our results. 
We can write the a gr as 



a 



9r -- + of + a 9r 



fini 



(4) 



where a 9r , a 9r and a 9 L n are the contributions in the soft, collinear and non-collinear regions. 
In the soft region, we can write the af as p2 



af = f dZ '{it) 

z ° V / bg 



of = a° 



a s r(l-e) Untf: 
2irT(l -2e) \ s 



(5) 



The lengthy expressions of the coefficients Aj will be given in Ref. [IB]. In the collinear 
region, the contributions after factorization can be written as two parts [22 



a 



■9r 



\\dzj bg 



'dL\ 9r fdL 
K dz) \dz 



gr 



(6) 



where the definition of the luminosity is similar to that in Eq. ||. Here the g and b are 



g/b(x,fi f ) = 



—fg/fb{x/y,fi f )Pgg /bb {y) 



(7) 



with 



where s is the subprocess center-of-mass energy and 



(8) 



1 + ? 2 

P bb (z) = C FJ ^ 
PJz) = -CM-z) 



P gg (z) = 2N 
P' gg {z) = 



z 1 — z ,„ . 
+ + z(l - z) 



l-z 



(9) 



with N = 3 and Cp = 4/3. & 9r can be written as 



a 



a s r(l - e) / Atx^I' 
2^T(1 - 2e) V s 



where 



-4f(fr ~+ frg) + jf(g ^ 99) + A , c(& 



bg)+A?{g^gg) , (10) 



4 SC 



Af(&^&<?) = Cir(21og5 s + 3/2) 

- <?<7) = 2N\og5 s + (UN - 2n,)/6 



(11) 
(12) 
(13) 



with ?2j = 5. 

For the initial-gluon processes, the results are much simpler compared to the gluon 
radiation processes, 



n ig — n ig _|_ n w 



(14) 



where cr* 9 and aj in are the contributions in the collinear and the non-collinear regions. After 



factorization, we can write al 9 as 22 



a 



dza c 



'dL\ 



ig _ 

z ° \ aZ ) gb 



(15) 



where the definition of b in the luminosity (jj^j - is 



b (x,fif)= —f g {x/y,iif)Pbg(y). 
Jx y 

The splitting functions in P bg , defined in Eq. [|, contains the parts 



(16) 



z 2 + a-z) 2 



Pl g {z) = -z{\ - z) 



(17) 



Another important issue, which does not exist in gluon-radiation processes, is the pro- 
cedure how to subtract the contribution of the on-shell anti-top quark decay to W~ and b 
from a t9 , besides subtracting double counting of g —>■ bb in collinear region. As in Ref. ||16|| , 
in order to remove all of the i contribution, we should subtract the term given by (in the 
narrow decay width limit) 



a = a LO {gg -> tt)B(t -> W'b), 



(18) 



where a LO (gg -> tt) and B(t -> W - 6) are LO cross section of (7(7 — > and branching ratio 
of the decay t — > W~b. 

It should be noticed that at 0(af), there are another QCD corrections arising from 
qq — > tbW~ and bq(q) — > [g = w, d, s], which can be treated by the similar methods 

described above. However, due to the lower luminosity of qq and bq(q) compared to those 
of gg and bg at the LHC, the QCD corrections are smaller, which are only a few percents. 
For completeness, we will include the initiated light quark contributions in our numerical 
results. 

Our numerical results are obtained using CTEQ5M (CTEQ5L) PDF |23| for NLO (LO) 
cross-section calculations. The 2-loop (1-loop) evolution of a s (fi) is adopted for NLO (LO) 
calculation and A^ 5 ^ = 226 (146) MeV for two-loop (one-loop) evolution. The top-quark 
pole mass is taken to be m t = 175 GeV; for simplicity, the bottom-quark mass has been 
omitted, and the renormalization and factorization scales are taken to be the same. We 



have compared the numerical results of the initial-gluon contribution to that in Ref. [IB 
and both results are in good agreement. 

In Fig. ^| we show the tree level and NLO cross sections as a function of renormalization 
and factorization scales /i//io (/-to = m t + mw). From the figure we can see that the NLO 
result is greater than the lowest order one. The NLO cross section is ~ 37 pb when \x = /io, 
while ~ 25 pb at tree level. In Fig. ^, the K factor (defined as the ratio of the NLO cross 
section to the LO one) is shown. For ^ < fi < 2//o, the K factor varies roughly between 
1.33 and 1.66. 

To summarize, the next-to-leading order QCD corrections to the Standard Model pro- 
cess bg — > tW~ at CERN large hadron collider with s/s = 14 TeV have been calculated. 
The NLO QCD corrections can enhance the cross section by a factor from 1.33 to 1.66 for 
renormalization and factorization scales ^ < /i < 2/jq. We should note here that the 
results presented in this letter are for the process bg — > tW~; they are the same for the 
charge conjugate process bg — > W + i. 
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FIG. 1. Feynmann diagrams for bg — ► : the Born level (a), the virtual gluon exchange (b) 
and the gluon radiation (c). 
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FIG. 2. Cross sections of NLO (dashed) and LO (solid) for pp — ► bg — > tW as functions of 
n/Ho at the LHC with yfs = 14 TeV, where = m t + my/- 
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FIG. 3. K factor for pp — > bg — ► tW as functions of /x//io at the LHC. 
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